Tox induces T cell IL-10 production in a BATF-dependent manner.

Tox is a member of the high mobility group (HMG)-Box transcription factors and plays important roles in thymic T cell development. Outside of the thymus, however, Tox is also highly expressed by CD8 and CD4 T cells in various states of activation and in settings of cancer and autoimmune disease. In CD4 T cells, Tox has been primarily studied in T follicular helper (TFH) cells where it, along with Tox2, promotes TFH differentiation by regulating key TFH-associated genes and suppressing CD4 cytotoxic T cell differentiation. However, the role of Tox in other T helper (Th) cell subtypes is less clear. Here, we show that Tox is expressed in several physiologically-activated Th subtypes and its ectopic expression enhances the in vitro differentiation of Th2 and T regulatory (Treg) cells. Tox overexpression in unpolarized Th cells also induced the expression of several genes involved in cell activation (Pdcd1), cellular trafficking (Ccl3, Ccl4, Xcl1) and suppressing inflammation (Il10) across multiple Th subtypes. We found that Tox binds the regulatory regions of these genes along with the transcription factors BATF, IRF4, and JunB and that Tox-induced expression of IL-10, but not PD-1, is BATF-dependent. Based on these data, we propose a model where Tox regulates Th cell chemotactic genes involved in facilitating dendritic cell-T cell interactions and aids in the resolution or prevention of inflammation through the production of IL-10.

Loss of CD4+ T cell-intrinsic arginase 1 accelerates Th1 response kinetics and reduces lung pathology during influenza infection.

Arginase 1 (Arg1), the enzyme catalyzing the conversion of arginine to ornithine, is a hallmark of IL-10-producing immunoregulatory M2 macrophages. However, its expression in T cells is disputed. Here, we demonstrate that induction of Arg1 expression is a key feature of lung CD4+ T cells during mouse in vivo influenza infection. Conditional ablation of Arg1 in CD4+ T cells accelerated both virus-specific T helper 1 (Th1) effector responses and its resolution, resulting in efficient viral clearance and reduced lung pathology. Using unbiased transcriptomics and metabolomics, we found that Arg1-deficiency was distinct from Arg2-deficiency and caused altered glutamine metabolism. Rebalancing this perturbed glutamine flux normalized the cellular Th1 response. CD4+ T cells from rare ARG1-deficient patients or CRISPR-Cas9-mediated ARG1-deletion in healthy donor cells phenocopied the murine cellular phenotype. Collectively, CD4+ T cell-intrinsic Arg1 functions as an unexpected rheostat regulating the kinetics of the mammalian Th1 lifecycle with implications for Th1-associated tissue pathologies.

C5a-licensed phagocytes drive sterilizing immunity during systemic fungal infection.

Systemic candidiasis is a common, high-mortality, nosocomial fungal infection. Unexpectedly, it has emerged as a complication of anti-complement C5-targeted monoclonal antibody treatment, indicating a critical niche for C5 in antifungal immunity. We identified transcription of complement system genes as the top biological pathway induced in candidemic patients and as predictive of candidemia. Mechanistically, C5a-C5aR1 promoted fungal clearance and host survival in a mouse model of systemic candidiasis by stimulating phagocyte effector function and ERK- and AKT-dependent survival in infected tissues. C5ar1 ablation rewired macrophage metabolism downstream of mTOR, promoting their apoptosis and enhancing mortality through kidney injury. Besides hepatocyte-derived C5, local C5 produced intrinsically by phagocytes provided a key substrate for antifungal protection. Lower serum C5a concentrations or a C5 polymorphism that decreases leukocyte C5 expression correlated independently with poor patient outcomes. Thus, local, phagocyte-derived C5 production licenses phagocyte antimicrobial function and confers innate protection during systemic fungal infection.

A comprehensive single cell data analysis of lymphoblastoid cells reveals the role of super-enhancers in maintaining EBV latency.

We probed the lifecycle of Epstein-Barr virus (EBV) on a cell-by-cell basis using single cell RNA sequencing (scRNA-seq) data from nine publicly available lymphoblastoid cell lines (LCLs). While the majority of LCLs comprised cells containing EBV in the latent phase, two other clusters of cells were clearly evident and were distinguished by distinct expression of host and viral genes. Notably, both were high expressors of EBV LMP1/BNLF2 and BZLF1 compared to another cluster that expressed neither gene. The two novel clusters differed from each other in their expression of EBV lytic genes, including glycoprotein gene GP350. The first cluster, comprising GP350- LMP1hi cells, expressed high levels of HIF1A and was transcriptionally regulated by HIF1-α. Treatment of LCLs with Pevonedistat, a drug that enhances HIF1-α signaling, markedly induced this cluster. The second cluster, containing GP350+ LMP1hi cells, expressed EBV lytic genes. Host genes that are controlled by super-enhancers (SEs), such as transcription factors MYC and IRF4, had the lowest expression in this cluster. Functionally, the expression of genes regulated by MYC and IRF4 in GP350+ LMP1hi cells were lower compared to other cells. Indeed, induction of EBV lytic reactivation in EBV+ AKATA reduced the expression of these SE-regulated genes. Furthermore, CRISPR-mediated perturbation of the MYC or IRF4 SEs in LCLs induced the lytic EBV gene expression, suggesting that host SEs and/or SE target genes are required for maintenance of EBV latency. Collectively, our study revealed EBV-associated heterogeneity among LCLs that may have functional consequence on host and viral biology.

Folic acid-mediated fibrosis is driven by C5a receptor 1-mediated activation of kidney myeloid cells.

We have previously reported that increased expression and activation of kidney cell complement components play an important role in the pathogenesis of renal scarring. Here, we used floxed green fluorescent protein (GFP)-C5a receptor 1 (C5aR1) knockin mice (GFP-C5ar1fl/fl) and the model of folic acid (FA)-induced kidney injury to define the cell types and potential mechanisms by which increased C5aR1 activation leads to fibrosis. Using flow cytometry and confocal microscopy, we identified macrophages as the major interstitial cell type showing increased expression of C5aR1 in FA-treated mice. C5ar1fl/fl.Lyz2Cre+/- mice, in which C5aR1 has been specifically deleted in lysozyme M-expressing myeloid cells, experienced reduced fibrosis compared with control C5ar1fl/fl mice. Examination of C5aR1-expressing macrophage transcriptomes by gene set enrichment analysis demonstrated that these cells were enriched in pathways corresponding to the complement cascade, collagen formation, and the NABA matrisome, strongly pointing to their critical roles in tissue repair/scarring. Since C5aR1 was also detected in a small population of platelet-derived growth factor receptor-ß+ GFP+ cells, we developed C5ar1fl/fl.Foxd1Cre+/- mice, in which C5aR1 is deleted specifically in pericytes, and found reduced FA-induced fibrosis. Primary cell cultures of platelet-derived growth factor receptor-ß+ pericytes isolated from FA-treated C5ar1fl/fl.Foxd1Cre+/- mice showed reduced secretion of several cytokines, including IL-6 and macrophage inflammatory protein-2, compared with pericytes isolated from FA-treated control GFP-C5ar1fl/fl mice. Collectively, these data imply that C5a/C5aR1 axis activation primarily in interstitial cells contributes to the development of renal fibrosis.NEW & NOTEWORTHY This study used novel green fluorescent protein C5a receptor 1 floxed mice and the model of folic acid-mediated kidney fibrosis to demonstrate the pathogenic role of increased expression of this complement receptor on macrophages.

Autocrine vitamin D signaling switches off pro-inflammatory programs of TH1 cells.

The molecular mechanisms governing orderly shutdown and retraction of CD4+ type 1 helper T (TH1) cell responses remain poorly understood. Here we show that complement triggers contraction of TH1 responses by inducing intrinsic expression of the vitamin D (VitD) receptor and the VitD-activating enzyme CYP27B1, permitting T cells to both activate and respond to VitD. VitD then initiated the transition from pro-inflammatory interferon-γ+ TH1 cells to suppressive interleukin-10+ cells. This process was primed by dynamic changes in the epigenetic landscape of CD4+ T cells, generating super-enhancers and recruiting several transcription factors, notably c-JUN, STAT3 and BACH2, which together with VitD receptor shaped the transcriptional response to VitD. Accordingly, VitD did not induce interleukin-10 expression in cells with dysfunctional BACH2 or STAT3. Bronchoalveolar lavage fluid CD4+ T cells of patients with COVID-19 were TH1-skewed and showed de-repression of genes downregulated by VitD, from either lack of substrate (VitD deficiency) and/or abnormal regulation of this system.

Host-Virus Chimeric Events in SARS-CoV-2-Infected Cells Are Infrequent and Artifactual.

The pathogenic mechanisms underlying severe SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) infection remain largely unelucidated. High-throughput sequencing technologies that capture genome and transcriptome information are key approaches to gain detailed mechanistic insights from infected cells. These techniques readily detect both pathogen- and host-derived sequences, providing a means of studying host-pathogen interactions. Recent studies have reported the presence of host-virus chimeric (HVC) RNA in transcriptome sequencing (RNA-seq) data from SARS-CoV-2-infected cells and interpreted these findings as evidence of viral integration in the human genome as a potential pathogenic mechanism. Since SARS-CoV-2 is a positive-sense RNA virus that replicates in the cytoplasm, it does not have a nuclear phase in its life cycle. Thus, it is biologically unlikely to be in a location where splicing events could result in genome integration. Therefore, we investigated the biological authenticity of HVC events. In contrast to true biological events like mRNA splicing and genome rearrangement events, which generate reproducible chimeric sequencing fragments across different biological isolates, we found that HVC events across >100 RNA-seq libraries from patients with coronavirus disease 2019 (COVID-19) and infected cell lines were highly irreproducible. RNA-seq library preparation is inherently error prone due to random template switching during reverse transcription of RNA to cDNA. By counting chimeric events observed when constructing an RNA-seq library from human RNA and spiked-in RNA from an unrelated species, such as the fruit fly, we estimated that ∼1% of RNA-seq reads are artifactually chimeric. In SARS-CoV-2 RNA-seq, we found that the frequency of HVC events was, in fact, not greater than this background "noise." Finally, we developed a novel experimental approach to enrich SARS-CoV-2 sequences from bulk RNA of infected cells. This method enriched viral sequences but did not enrich HVC events, suggesting that the majority of HVC events are, in all likelihood, artifacts of library construction. In conclusion, our findings indicate that HVC events observed in RNA-sequencing libraries from SARS-CoV-2-infected cells are extremely rare and are likely artifacts arising from random template switching of reverse transcriptase and/or sequence alignment errors. Therefore, the observed HVC events do not support SARS-CoV-2 fusion to cellular genes and/or integration into human genomes. IMPORTANCE The pathogenic mechanisms underlying SARS-CoV-2, the virus responsible for COVID-19, are not fully understood. In particular, relatively little is known about the reasons some individuals develop life-threatening or persistent COVID-19. Recent studies identified host-virus chimeric (HVC) reads in RNA-sequencing data from SARS-CoV-2-infected cells and suggested that HVC events support potential "human genome invasion" and "integration" by SARS-CoV-2. This suggestion has fueled concerns about the long-term effects of current mRNA vaccines that incorporate elements of the viral genome. SARS-CoV-2 is a positive-sense, single-stranded RNA virus that does not encode a reverse transcriptase and does not include a nuclear phase in its life cycle, so some doubts have rightfully been expressed regarding the authenticity of HVCs and the role played by endogenous retrotransposons in this phenomenon. Thus, it is important to independently authenticate these HVC events. Here, we provide several lines of evidence suggesting that the observed HVC events are likely artifactual.

SARS-CoV-2 drives JAK1/2-dependent local complement hyperactivation.

Patients with coronavirus disease 2019 (COVID-19) present a wide range of acute clinical manifestations affecting the lungs, liver, kidneys and gut. Angiotensin converting enzyme (ACE) 2, the best-characterized entry receptor for the disease-causing virus SARS-CoV-2, is highly expressed in the aforementioned tissues. However, the pathways that underlie the disease are still poorly understood. Here, we unexpectedly found that the complement system was one of the intracellular pathways most highly induced by SARS-CoV-2 infection in lung epithelial cells. Infection of respiratory epithelial cells with SARS-CoV-2 generated activated complement component C3a and could be blocked by a cell-permeable inhibitor of complement factor B (CFBi), indicating the presence of an inducible cell-intrinsic C3 convertase in respiratory epithelial cells. Within cells of the bronchoalveolar lavage of patients, distinct signatures of complement activation in myeloid, lymphoid and epithelial cells tracked with disease severity. Genes induced by SARS-CoV-2 and the drugs that could normalize these genes both implicated the interferon-JAK1/2-STAT1 signaling system and NF-κB as the main drivers of their expression. Ruxolitinib, a JAK1/2 inhibitor, normalized interferon signature genes and all complement gene transcripts induced by SARS-CoV-2 in lung epithelial cell lines, but did not affect NF-κB-regulated genes. Ruxolitinib, alone or in combination with the antiviral remdesivir, inhibited C3a protein produced by infected cells. Together, we postulate that combination therapy with JAK inhibitors and drugs that normalize NF-κB-signaling could potentially have clinical application for severe COVID-19.

Host-virus chimeric events in SARS-CoV2 infected cells are infrequent and artifactual.

Pathogenic mechanisms underlying severe SARS-CoV2 infection remain largely unelucidated. High throughput sequencing technologies that capture genome and transcriptome information are key approaches to gain detailed mechanistic insights from infected cells. These techniques readily detect both pathogen and host-derived sequences, providing a means of studying host-pathogen interactions. Recent studies have reported the presence of host-virus chimeric (HVC) RNA in RNA-seq data from SARS-CoV2 infected cells and interpreted these findings as evidence of viral integration in the human genome as a potential pathogenic mechanism. Since SARS-CoV2 is a positive sense RNA virus that replicates in the cytoplasm it does not have a nuclear phase in its life cycle, it is biologically unlikely to be in a location where splicing events could result in genome integration. Here, we investigated the biological authenticity of HVC events. In contrast to true biological events such as mRNA splicing and genome rearrangement events, which generate reproducible chimeric sequencing fragments across different biological isolates, we found that HVC events across >100 RNA-seq libraries from patients with COVID-19 and infected cell lines, were highly irreproducible. RNA-seq library preparation is inherently error-prone due to random template switching during reverse transcription of RNA to cDNA. By counting chimeric events observed when constructing an RNA-seq library from human RNA and spike-in RNA from an unrelated species, such as fruit-fly, we estimated that ~1% of RNA-seq reads are artifactually chimeric. In SARS-CoV2 RNA-seq we found that the frequency of HVC events was, in fact, not greater than this background "noise". Finally, we developed a novel experimental approach to enrich SARS-CoV2 sequences from bulk RNA of infected cells. This method enriched viral sequences but did not enrich for HVC events, suggesting that the majority of HVC events are, in all likelihood, artifacts of library construction. In conclusion, our findings indicate that HVC events observed in RNA-sequencing libraries from SARS-CoV2 infected cells are extremely rare and are likely artifacts arising from either random template switching of reverse-transcriptase and/or sequence alignment errors. Therefore, the observed HVC events do not support SARS-CoV2 fusion to cellular genes and/or integration into human genomes.

CDK8 Fine-Tunes IL-6 Transcriptional Activities by Limiting STAT3 Resident Time at the Gene Loci.

Cytokines are highly pleiotropic ligands that regulate the immune response. Here, using interleukin-6 (IL-6) as a model system, we perform detailed phosphoproteomic and transcriptomic studies in human CD4+ T helper 1 (Th-1) cells to address the molecular bases defining cytokine functional pleiotropy. We identify CDK8 as a negative regulator of STAT3 transcriptional activities, which interacts with STAT3 upon IL-6 stimulation. Inhibition of CDK8 activity, using specific small molecule inhibitors, reduces the IL-6-induced phosphoproteome by 23% in Th-1 cells, including STAT3 S727 phosphorylation. STAT3 binding to target DNA sites in the genome is increased upon CDK8 inhibition, which results in a concomitant increase in STAT3-mediated transcriptional activity. Importantly, inhibition of CDK8 activity under Th-17 polarizing conditions results in an enhancement of Th-17 differentiation. Our results support a model where CDK8 regulates STAT3 transcriptional processivity by modulation of its gene loci resident time, critically contributing to diversification of IL-6 responses.

Epstein-Barr Virus Episome Physically Interacts with Active Regions of the Host Genome in Lymphoblastoid Cells.

The Epstein-Barr virus (EBV) episome is known to interact with the three-dimensional structure of the human genome in infected cells. However, the exact locations of these interactions and their potential functional consequences remain unclear. Recently, high-resolution chromatin conformation capture (Hi-C) assays in lymphoblastoid cells have become available, enabling us to precisely map the contacts between the EBV episome(s) and the human host genome. Using available Hi-C data at a 10-kb resolution, we have identified 15,000 reproducible contacts between EBV episome(s) and the human genome. These contacts are highly enriched in chromatin regions denoted by typical or super enhancers and active markers, including histone H3K27ac and H3K4me1. Additionally, these contacts are highly enriched at loci bound by host transcription factors that regulate B cell growth (e.g., IKZF1 and RUNX3), factors that enhance cell proliferation (e.g., HDGF), or factors that promote viral replication (e.g., NBS1 and NFIC). EBV contacts show nearly 2-fold enrichment in host regions bound by EBV nuclear antigen 2 (EBNA2) and EBNA3 transcription factors. Circular chromosome conformation capture followed by sequencing (4C-seq) using the EBV origin of plasmid replication (oriP) as a "bait" in lymphoblastoid cells further confirmed contacts with active chromatin regions. Collectively, our analysis supports interactions between EBV episome(s) and active regions of the human genome in lymphoblastoid cells.IMPORTANCE EBV is associated with ∼200,000 cancers each year. In vitro, EBV can transform primary human B lymphocytes into immortalized cell lines. EBV-encoded proteins, along with noncoding RNAs and microRNAs, hijack cellular proteins and pathways to control cell growth. EBV nuclear proteins usurp normal transcriptional programs to activate the expression of key oncogenes, including MYC, to provide a proliferation signal. EBV nuclear antigens also repress CDKN2A to suppress senescence. EBV membrane protein activates NF-κB to provide survival signals. EBV genomes are maintained by EBNA1, which tethers EBV episomes to the host chromosomes during mitosis. However, little is known about where EBV episomes are located in interphase cells. In interphase cells, EBV promoters drive the expression of latency genes, while oriP functions as an enhancer for these promoters. In this study, integrative analyses of published lymphoblastoid cell line (LCL) Hi-C data and our 4C-seq experiments position EBV episomes to host genomes with active epigenetic marks. These contact points were significantly enriched for super enhancers. The close proximity of EBV episomes and the super enhancers that are enriched for transcription cofactors or mediators in lymphoblasts may benefit EBV gene expression, suggesting a novel mechanism of transcriptional activation.

An autocrine Vitamin D-driven Th1 shutdown program can be exploited for COVID-19.

Pro-inflammatory immune responses are necessary for effective pathogen clearance, but cause severe tissue damage if not shut down in a timely manner 1,2 . Excessive complement and IFN-γ-associated responses are known drivers of immunopathogenesis 3 and are among the most highly induced immune programs in hyper-inflammatory SARS-CoV2 lung infection 4 . The molecular mechanisms that govern orderly shutdown and retraction of these responses remain poorly understood. Here, we show that complement triggers contraction of IFN-γ producing CD4 + T helper (Th) 1 cell responses by inducing expression of the vitamin D (VitD) receptor (VDR) and CYP27B1, the enzyme that activates VitD, permitting T cells to both activate and respond to VitD. VitD then initiates the transition from pro-inflammatory IFN-γ + Th1 cells to suppressive IL-10 + Th1 cells. This process is primed by dynamic changes in the epigenetic landscape of CD4 + T cells, generating superenhancers and recruiting c-JUN and BACH2, a key immunoregulatory transcription factor 5-7 . Accordingly, cells in psoriatic skin treated with VitD increased BACH2 expression, and BACH2 haplo-insufficient CD4 + T cells were defective in IL-10 production. As proof-of-concept, we show that CD4 + T cells in the bronchoalveolar lavage fluid (BALF) of patients with COVID-19 are Th1-skewed and that VDR is among the top regulators of genes induced by SARS-CoV2. Importantly, genes normally down-regulated by VitD were de-repressed in CD4 + BALF T cells of COVID-19, indicating that the VitD-driven shutdown program is impaired in this setting. The active metabolite of VitD, alfacalcidol, and cortico-steroids were among the top predicted pharmaceuticals that could normalize SARS-CoV2 induced genes. These data indicate that adjunct therapy with VitD in the context of other immunomodulatory drugs may be a beneficial strategy to dampen hyperinflammation in severe COVID-19.

SARS-CoV2 drives JAK1/2-dependent local and systemic complement hyper-activation.

Patients with coronavirus disease 2019 (COVID-19) present with a range of devastating acute clinical manifestations affecting the lungs, liver, kidneys and gut. The best-characterized entry receptor for the disease-causing virus SARS-CoV2, angiotensin converting enzyme (ACE) 2, is highly expressed in these tissues. However, the pathways that underlie the disease are still poorly understood. Here we show that the complement system is unexpectedly one of the intracellular pathways most highly induced by SARS-CoV2 infection in lung epithelial and liver cells. Within cells of the bronchoalveolar lavage of patients, distinct signatures of complement activation in myeloid, lymphoid and epithelial cells tracked with disease severity. Modelling the regulome of host genes induced by COVID-19 and the drugs that could normalize these genes both implicated the JAK1/2-STAT1 signaling system downstream of type I interferon receptors, and NF-kB. Ruxolitinib, a JAK1/2 inhibitor and the top predicted pharmaceutical candidate, normalized interferon signature genes, IL-6 (the best characterized severity marker in COVID-19) and all complement genes induced by SARS-CoV2, but did not affect NF-kB-regulated genes. We predict that combination therapy with JAK inhibitors and other agents with the potential to normalize NF-kB-signaling, such as anti-viral agents, may serve as an effective clinical strategy.

Diapedesis-Induced Integrin Signaling via LFA-1 Facilitates Tissue Immunity by Inducing Intrinsic Complement C3 Expression in Immune Cells.

Intrinsic complement C3 activity is integral to human T helper type 1 (Th1) and cytotoxic T cell responses. Increased or decreased intracellular C3 results in autoimmunity and infections, respectively. The mechanisms regulating intracellular C3 expression remain undefined. We identified complement, including C3, as among the most significantly enriched biological pathway in tissue-occupying cells. We generated C3-reporter mice and confirmed that C3 expression was a defining feature of tissue-immune cells, including T cells and monocytes, occurred during transendothelial diapedesis, and depended on integrin lymphocyte-function-associated antigen 1 (LFA-1) signals. Immune cells from patients with leukocyte adhesion deficiency type 1 (LAD-1) had reduced C3 transcripts and diminished effector activities, which could be rescued proportionally by intracellular C3 provision. Conversely, increased C3 expression by T cells from arthritis patients correlated with disease severity. Our study defines integrins as key controllers of intracellular complement, demonstrates that perturbations in the LFA-1-C3-axis contribute to primary immunodeficiency, and identifies intracellular C3 as biomarker of severity in autoimmunity.

Kinetics of cytokine receptor trafficking determine signaling and functional selectivity.

Cytokines activate signaling via assembly of cell surface receptors, but it is unclear whether modulation of cytokine-receptor binding parameters can modify biological outcomes. We have engineered IL-6 variants with different affinities to gp130 to investigate how cytokine receptor binding dwell-times influence functional selectivity. Engineered IL-6 variants showed a range of signaling amplitudes and induced biased signaling, with changes in receptor binding dwell-times affecting more profoundly STAT1 than STAT3 phosphorylation. We show that this differential signaling arises from defective translocation of ligand-gp130 complexes to the endosomal compartment and competitive STAT1/STAT3 binding to phospho-tyrosines in gp130, and results in unique patterns of STAT3 binding to chromatin. This leads to a graded gene expression response and differences in ex vivo differentiation of Th17, Th1 and Treg cells. These results provide a molecular understanding of signaling biased by cytokine receptors, and demonstrate that manipulation of signaling thresholds is a useful strategy to decouple cytokine functional pleiotropy.

Integrated Pan-Cancer Map of EBV-Associated Neoplasms Reveals Functional Host-Virus Interactions.

Epstein-Barr virus (EBV) is a complex oncogenic symbiont. The molecular mechanisms governing EBV carcinogenesis remain elusive and the functional interactions between virus and host cells are incompletely defined. Here we present a comprehensive map of the host cell-pathogen interactome in EBV-associated cancers. We systematically analyzed RNA sequencing from >1,000 patients with 15 different cancer types, comparing virus and host factors of EBV+ to EBV- tissues. EBV preferentially integrated at highly accessible regions of the cancer genome, with significant enrichment in super-enhancer architecture. Twelve EBV transcripts, including LMP1 and LMP2, correlated inversely with EBV reactivation signature. Overexpression of these genes significantly suppressed viral reactivation, consistent with a "virostatic" function. In cancer samples, hundreds of novel frequent missense and nonsense variations in virostatic genes were identified, and variant genes failed to regulate their viral and cellular targets in cancer. For example, one-third of patients with EBV+ NK/T-cell lymphoma carried two novel nonsense variants (Q322X, G342X) of LMP1 and both variant proteins failed to restrict viral reactivation, confirming loss of virostatic function. Host cell transcriptional changes in response to EBV infection classified tumors into two molecular subtypes based on patterns of IFN signature genes and immune checkpoint markers, such as PD-L1 and IDO1. Overall, these findings uncover novel points of interaction between a common oncovirus and the human genome and identify novel regulatory nodes and druggable targets for individualized EBV and cancer-specific therapies. SIGNIFICANCE: This study provides a comprehensive map of the host cell-pathogen interactome in EBV+ malignancies.See related commentary by Mbulaiteye and Prokunina-Olsson, p. 5917.

PRMDA: personalized recommendation-based MiRNA-disease association prediction.

Recently, researchers have been increasingly focusing on microRNAs (miRNAs) with accumulating evidence indicating that miRNAs serve as a vital role in various biological processes and dysfunctions of miRNAs are closely related with human complex diseases. Predicting potential associations between miRNAs and diseases is attached considerable significance in the domains of biology, medicine, and bioinformatics. In this study, we developed a computational model of Personalized Recommendation-based MiRNA-Disease Association prediction (PRMDA) to predict potential related miRNA for all diseases by implementing personalized recommendation-based algorithm based on integrated similarity for diseases and miRNAs. PRMDA is a global method capable of prioritizing candidate miRNAs for all diseases simultaneously. Moreover, the model could be applied to diseases without any known associated miRNAs. PRMDA obtained AUC of 0.8315 based on leave-one-out cross validation, which demonstrated that PRMDA could be regarded as a reliable tool for miRNA-disease association prediction. Besides, we implemented PRMDA on the HMDD V1.0 and HMDD V2.0 databases for three kinds of case studies about five important human cancers in order to test the performance of the model from different perspectives. As a result, 92%, 94%, 88%, 96% and 88% out of the top 50 candidate miRNAs predicted by PRMDA for Colon Neoplasms, Esophageal Neoplasms, Lymphoma, Lung Neoplasms and Breast Neoplasms, respectively, were confirmed by experimental reports.